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Abstract 

We construct six multi-parameter families of Hermitian quasi-exactly solvable matrix 
Schrodinger operators in one variable. The method for finding these operators relies heavily 
upon a special representation of the Lie algebra o(2, 2) = sl(2) sl(2) whose representation 
CN ' space contains an invariant finite-dimensional subspace. Besides that we give several ex- 

amples of quasi-exactly solvable matrix models that have square-integrable eigenfunctions. 
These examples are in direct analogy with the quasi-exactly solvable scalar Schrodinger 



operators obtained by Turbiner and Ushveridze. 



I. Introduction 



In the paper we have suggested a generalization of the Turbiner-Shifman approach @-@| to 
the construction of quasi-exactly solvable (QES) models on line for the case of matrix Hamilto- 
^ ■ nians. We remind that originally their method was applied to scalar one-dimensional stationary 
Schrodinger equations. Later on it was extended to the case of mult i- dimensional scalar station- 
ary Schrodinger equations 0-0 (see, also ||). 

A systematic description of our approach can be found in the paper f9j. The procedure of 
constructing a QES matrix (scalar) model is based on the concept of a Lie-algebraic Hamiltonian. 
We call a second-order operator in one variable Lie-algebraic if the following requirements are 
met: 

• The Hamiltonian is a quadratic form with constant coefficients of first-order operators 
Qi, Q2, ■ ■ ■ ,Q n forming a Lie algebra g\ 

• The Lie algebra g has a finite-dimensional invariant subspace X of the whole representation 
space. 

Now if a given Hamiltonian H[x] is Lie-algebraic, then after being restricted to the space X 
it becomes a matrix operator TC whose eigenvalues and eigenvectors are computed in a purely 
algebraic way. This means that the Hamiltonian H[x] is quasi-exactly solvable (for further details 
on scalar QES models see 0). 
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It should be noted that there exist alternative approaches to constructing matrix QES models 
K|-(l5|]. The principal idea of these is fixing the form of basis elements of the invariant space 



X. They are chosen to be polynomials in x. This assumption leads to a challenging problem of 



classification of superalgebras by matrix-differential operators in one variable |L5 . 

We impose no a priori restrictions on the form of basis elements of the space I. What is 
fixed it is the class to which the basis elements of the Lie algebra g should belong. Following 
0, H we choose this class £ as the set of matrix differential operators of the form 

£ = {Q: Q = a(x)d x + A(x)} . (1) 

Here a(x) is a smooth real-valued function and A(x) is an iV x iV matrix whose entries are 
smooth complex- valued functions of x. Hereafter we denote d/dx as d x . 

Evidently, £ can be treated as an infinite-dimensional Lie algebra with a standard com- 
mutator as a Lie bracket. Given a subalgebra (Q±, Q2, ■ ■ ■ , Q n ) of the algebra £, that has a 
finite-dimensional invariant space, we can easily construct a QES matrix model. To this end we 
compose a bilinear combination of the operators Qi, Q2, ■ ■ ■ , Q n and of the unit N x N matrix 
/ with constant complex coefficients ctjk and get 

H[x]ip(x) = f jr a 3k Q 3 Q k \ . (2) 

So there arises a natural problem of classification of subalgebras of the algebra £ within 
its inner automorphism group. The problem of classification of inequivalent realizations of Lie 
algebras on line and on plane has been solved in a full generality by Lie itself [16|, [17] (see, also 



fl8|j). However, the classification problem for the case when A(x) 7^ f(x)I with a scalar function 
f(x) is open by now. In the paper || we have classified realizations of the Lie algebras of the 
dimension up to three by the operators belonging to £ with an arbitrary N. Next, fixing N = 2 
we have studied which of them give rise to QES matrix Hamiltonians H[x]. It occurs that the 
only three-dimensional algebra that meets this requirement is the algebra s/(2) (which is fairly 
easy to predict taking into account the scalar case!). This yields the two families of 2 x 2 QES 
models, one of them under proper restrictions giving rise to the well-known family of scalar QES 
Hamiltonians (for more details, see |9|]). 

As is well-known a physically meaningful QES matrix Schrodinger operator has to be Hermi- 
tian. This requirement imposes restrictions on the choice of QES models which somehow were 
beyond considerations of our previous papers Jl|, [|. The principal aim of the present paper is 
to formulate and implement an efficient algebraic procedure for constructing QES Hermitian 
matrix Schrodinger operators 

H[x] = d 2 x + V(x). (3) 

This requires a slight modification of the algebraic procedure used in ||. We consider as an 
algebra g the direct sum of two s/(2) algebras which is equivalent to the algebra o(2,2). The 
necessary algebraic structures are introduced in Section 2. The next Section is devoted to 
constructing in a regular way Hermitian QES matrix Schrodinger operators on line. We give the 
list of thus obtained QES models in Section 4. The fifth Section contains a number of examples 
of Hermitian QES Schrodinger operators that have square integrable eigenfunctions. 
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(4) 



II. Extension of the algebra si (2) 

Following |], [| we consider the representation of the algebra sl(2) 

sl(2) = (Q_, Q , Q+) 

= (d x , xd x - Z±f± + S , x 2 d x - (m - l)x + 2S x + S + ), 

where So = cr 3 /2, S + = [io 2 + o"i)/2, o~fc are the 2x2 Pauli matrices 

ai = ( °i I ) ' ° 2 = ( ? o* ) ' as = ( I -1 ) 

and m > 2 is an arbitrary natural number. This representation gives rise to a family of QES 
models and furthermore the algebra (f|) has the following finite-dimensional invariant space 

1 s i{2) = Xi 0X 2 = (ei, xei, . . . , x m ~ 2 e x ) 

(5) 

(me*2, . . . , mx^e.2 — jx^~ e*i, . . . , mx m e*2 — mx m_1 e*i). 

Since the spaces Xi, X2 are invariant with respect to an action of any of the operators (^), the 
above representation is reducible. A more serious trouble is that it is not possible to construct 
a QES operator, that is equivalent to a Hermitian Schrodinger operator, by taking a bilinear 
combination (|2|) of operators (|]) with coefficients being complex numbers. To overcome this 
difficulty we use the idea indicated in || and let the coefficients of the bilinear combination @j 
to be constant 2x2 matrices. To this end we introduce a wider Lie algebra and add to the 
algebra @) the following three matrix operators: 

R- = S-, Ro = S-X + Sq, R+ = S^x 2 + 2Sqx + S+, (6) 

where S± = (io<i ± o\)j2. 

It is straightforward to verify that the space (||) is invariant with respect to an action of a 
linear combination of the operators (0). Consider next the following set of operators: 

(T± = Q± — R±, T = Q — R , R±, R , I), (7) 

where Q and R are operators (^) and (Q), respectively, and J is a unit 2x2 matrix. By a 
direct computation we check that the operators T±, T as well as the operators R±, Ro, fulfill the 
commutation relations of the algebra sl(2). Furthermore any of the operators T±,Tq commutes 
with any of the operators R±,Rq. Consequently, operators (0) form the Lie algebra 

S /(2)0 S /(2) 0/-o(2,2)0J. 

In a sequel we denote this algebra as g. 

The Casimir operators of the Lie algebra g are multiples of the unit matrix 

d = T 2 - T + T_ - To = [-^) 1, K 2 = R 2 Q - R + R_ - R = -I. 
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Using this fact it can be shown that the representation of g realized on the space 2 s i(2) is 
irreducible. 

One more remark is that the operators (0) satisfy the following relations: 
R 2 _ = 0, Rq = 2j;, R 2 + = 0, 

{R^ 1 R } = 1 {R +1 R } = 1 {R„,R + } = -l, (8) 
R—Rq = ijR—, RqR+ = 2^?+) R—R+ = Rq — 2- 

Here {Qi, Q2} = Q1Q2 + Q2Q\- One of the consequences of this fact is that the algebra g may 
be considered as a superalgebra which shows and evident link to the results of the paper JT5 . 



III. The general form of the Hermitian QES operator 

Using the commutation relations of the Lie algebra g together with relations (^) one can show 
that any bilinear combination of the operators (|7|) is a linear combination of twenty one (basis) 
quadratic forms of the operators (0). Composing this linear combination yields all QES models 
which can be obtained with the help of our approach. However the final goal of the paper is not 
to get some families of QES matrix second-order operators as such but to get QES Schrodinger 
operators ([3]). This means that it is necessary to transform bilinear combination (^) to the 
standard form ([3]). What is more, it is essential that the corresponding transformation should 
be given by explicit formulae, since we need to write down explicitly the matrix potential V(x) 
of thus obtained QES Schrodinger operator and the basis functions of its invariant space. 

The general form of QES model obtainable within the framework of our approach is as follows 

H[x] = i(x)d 2 x + B(x)d x + C(x), (9) 

where £(x) is some real- valued function and B(x),C(x) are matrix functions of the dimension 
2x2. Let U(x) be an invertible 2x2 matrix-function satisfying the system of ordinary differential 
equations 

U'(x) = ^ - B(x)^j U(x), (10) 

and the function f(x) be defined by the relation 



Then the change of variables reducing @ to the standard form (Q) reads as 

X - „ = /(«), 

H[x] -. H[y\ = U-\y)H[r\y)\U(y), 

where f~ l stands for the inverse of / and U(y) = U(f~ l (y)). 

Performing the transformation (|T2"| ) yields the Schrodinger operator 

H[y] = dl + V(y) (13) 
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with 



V(y) 



u-h 



x 



±B 2 (x) - \B'{x) + ^B(x) + C(x) 



U(x) 



e _ 3^ 



l2 



16? 



(14) 



x=f-Hy) 

Hereafter, the notation {W(x)} x= f-i( y ) means that we should replace x with f~ 1 (y) in the 
expression W(x). 

Furthermore, if we denote the basis elements of the invariant space @ as fi(x), . . . , / 2m (i), 
then the invariant space of the operator H[y] takes the form 

(15) 



lsi (2) = (U- 1 (y)f 1 (f- 1 (y)) 1 . . . , U-\y)f 2m {r\y))) ■ 



In view of the remark made at the beginning of this section we are looking for such QES 
models that the transformation law ([12]) can be given explicitly. This means that we should be 
able to construct a solution of system ([It]) in an explicit form. To achieve this goal we select 
from the above mentioned set of twenty one linearly independent quadratic forms of operators 
(0) the twelve forms, 



A 



A l = xd 2 , A 2 = x 2 d 2 x + (m - l)a 3 
°3 „ „ 2 



d x , B 1 = xd x + 



Bo 



x d x — {m — l)x + a 3 x + a x , 



a , m 

VlO x + y<7 3 , 



m 



<?zd x 



(16) 



C 2 = ia 2 d x + —cr 3 , C 3 

B x = x 3 dl - 2a x xd x + (3m - m 2 - 3)x + (2m - 3)xa 3 + (4m - A)a x , 
D 2 = x 3 d^ — 2ia 2 xd x + (3m — m 2 — 3)x + (2m — 3)xa 3 + (4m — 4)<7i, 
D 3 = 2a 3 xd x + (1 - 2m) a 3 , 

whose linear combinations have such a structure that system (|T0|) , can be integrated in a closed 
form. However, in the present paper we study systematically the first nine quadratic forms from 
the above list. The quadratic forms D X ,D 2 , D 3 are used to construct an example of QES model 
such that the matrix potential is expressed via the Weierstrass function. 

Thus the general form of the Hamiltonian to be considered in a sequel is as follows 



H\x\ 



2 3 

^(a M A M + Pf.B^) + ^7iCi = (a 2 x 2 + a x x + a Q )d 2 x 

/i=0 i=l 

+ (P 2 x 2 + (3 x x + A) + TiO"! + ij 2 cr 2 + 730-3)^ + (3 2 a 3 x 



(17) 



—/3 2 (m - l)x + (3 2 oi + 



a 2 [m 



1) + y + y(7i + 7 2 ) 



o"3- 



Here a , a x , a 2 are arbitrary real constants and (3 , 
If we denote 



, 73 are arbitrary complex constants. 



7i = 7i> 72 = 172, 73 = 73, $ = 2a 2 (m - 1) + f3 x + m(7i + 72) 
£(x) = a 2 x 2 + a\x + «o, v( x ) = P2X 2 + /3 x x + /3q, 
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then the general solution of system ( JlQ| ) reads as 



U(x) = £ 1/4 (x)exp 



Tj(x) 



dx 



cxp 



-dx 



A. 



(19) 



where A is an arbitrary constant invertible 2x2 matrix. Performing the transformation (ji 
with U(x) being given by fll~9l ) reduces QES operator (|I7D to a Schrodinger form (|I~3D , where 



/Ia"^ V + 2C't7 - 2£t/ - 4/3 2 (m - l)x£ - 7 - 
+2(f - 77)7itTi + ^U-\x)a x U{x) + (A(3 2 x + 25)£ 



(20) 



TT-if \ tti m a , a 2 3(2a 2 x + «i) 2 
XC/ i (x)a 3 f/(x)}A + y-^ — ^ 



Here £,17 are functions of x defined in fll8|) and / 1 (y) is the inverse of f(x) which is given by 

(0)- 

The requirement of hermiticity of the Schrodinger operator ( jl3|) is equivalent to the require- 
ment of hermiticity of the matrix V(y). To select from the multi-parameter family of matrices 
([ZD]) Hermitian ones we will make use of the following technical lemmas. 

Lemma 1 The matrices za a ,w(a a ± ia^),a ^ b, with {z, w} C C,z R, w ^ cannot be 
reduced to Hermitian matrices with the help of a transformation 



A^ A' = A~ l AA, 
where A is an invertible constant 2x2 matrix. 



(21) 



Proof. It is sufficient to prove the statement for the case a = 1, b = 2, since all other cases 



are equivalent to this one. Suppose the inverse, namely that there exists a transformation (pi ) 
transforming the matrix zo\ to a Hermitian matrix A'. As tr (zo~i) = tr A' = 0, the matrix A' 
has the form with some real constants a«. Next, from the equality det (zai) = det A' we get 
z 2 = a 2 . The last relation is in contradiction to the fact that z ^ R. Consequently, the matrix 
zo\ cannot be reduced to a Hermitian matrix with the aid of a transformation ([H]) . 
Let us turn now to the matrix w{a\ + io-i). Taking a general form of the matrix A 

A = 

we represent (|2T|) as follows 

A' = A' 1 w(a 1 + ia 2 )A = ^ ( ^ J , 5 = det A. 

The conditions of hermiticity of the matrix A' read 

W W - —W 2 W -2 

-red = ■=■ ca, — — c = ■=■ a . 
oo oo 
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where the bar over a symbol stands for the complex conjugation. 

It follows from the second relation that c, d can vanish only simultaneously which is impossible 
in view of the fact that the matrix A is invertible. Consequently, the relation cd 7^ holds. Hence 
we get 

Z^ = §^| C |2 + M |2 = . 

c a 

This contradiction proves the fact that the matrix w(a\ + icr 2 ) cannot be reduced to a Hermitian 
form. 

As the matrix <j\ + ioi is transformed to become <j\ — ia 2 with the use of an appropriate 
transformation (|2l|), the lemma is proved. 



Lemma 2 Let a = (01,02,03), b = (61,62,^3), c = (01,02,03) be complex vectors and a be the 
vector whose components are the Pauli matrices {a\, 0%, 03)- Then the following assertions holds 
true. 

1. A non-zero matrix aa is reduced to a Hermitian form with the help of a transformation 
iff' a 2 > ( this inequality means, in particular, that a 2 G R ); 



2. Non-zero matrices aa, ba with b 7^ Aa, A G R, are reduced simultaneously to Hermitian 
forms with the help of a transformation $2\ ) iff 



a 2 > 0, ? > 0, (ax b) 2 > 0; 

3. Matrices aa, ba, ca with a 7^ 0, 67^ Aa, fib, {A, //} C R are reduced simultaneously to 
Hermitian forms with the help of a transformation $Z1\ ) iff 

a 2 > 0, b 2 > 0, (a x b) 2 > 0, 
|ac, be, (a xb)c\ C R. 

Here we designate the scalar product of vectors a, b as ah and the vector product of these as axb. 

Proof. Let us first prove the necessity of the assertion 1 of the lemma. Suppose that the 
non-zero matrix aa can be reduced to a Hermitian form. We will prove that hence it follows the 
inequality a 2 > 0. 

Consider the matrices: 

, 1 + eiik^- ^J^ - h ia k , a^0, 
Kj{a,b) = { J a T ' (22) 

1, a = 0, 

where (i,j,k) = cycle (1,2,3). It is not difficult to verify that these matrices are invertible, 
provided 



a 2 + a} ^ 0. (23) 
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Given the condition the following relations hold 



A i - 1 (a,6) (X/ Aij(a,b) 



hoi + acTj 
Va 2 + 6 2 



6(7, 



Va 2 + b 2 



I = k. 



(24) 



As a is a non-zero vector, there exists at least one pair of the indices i, j such that af + a 2 ^ 0. 



Applying the transformation 



with a = cii,b = dj we get 



aa 



a a 



a 2 + a 2 Oj + a k a k 



(25) 



(no summation over the indices k is carried out). As the direct check shows, the quantity a 2 
is invariant with respect to transformation 



i.e. a 



?/2 



If a 2 = 0, then a' 2 + a!{ = 0, or a\ = ±ia' k . Hence by force of Lemma 1 it follows that 
the matrix (p5|) cannot be reduced to a Hermitian form. Consequently, a 2 ^ and the relation 

]) with a = 



a ' 2 j + 7^ holds true. Applying transformation 



a 2 + a 2 , b = a k we get 



a a 



(26) 



Due to Lemma 1, if the number a/^ 2 is complex, then the above matrix cannot be transformed 
to a Hermitian matrix. Consequently, the relation a 2 > holds true. 

The sufficiency of the assertion 1 of the lemma follows from the fact that, given the condition 
a 2 > 0, the matrix (|2B|) is Hermitian. 

Now we will prove the necessity of the assertion 2 of the lemma. First of all we note that 
due to assertion 1, a 2 >,b 2 > 0. Next, without loss of generality we can again suppose that 



a? 



+ a 2 ^ 0. Taking the superposition of two transformations of the form (|24| ) with a 



a,; , b 



and a = Jaf + aj,b = a k yields 



Kj(a-i, a-j)Ajk(^a 2 + a 2 , a k ) = 1 + ie ijk v , a 2 = a ^ 



a 2 



a; 



-<7 k t€ij k 



a 2 - ak c 
a 2 + a 2 



(27) 



(here the finite limit exists when — >• 0). Using this formula and taking into account (|24|) yield 



acr 



6a — > 6'c? 



6jOj 



a| + a 2 



6,-Oj a k ab — b k a 2 

-CTi + 



(P s /a 2 + a 2 3 



ab 



'a 2 



Ok- 



:28i 



Let us show that the necessary condition for the matrices y3?a kl b'a to be reducible to Her- 
mitian forms simultaneously reads as ab e R. Indeed, as the matrices b'a, a k are simultaneously 
reduced to Hermitian forms, the matrix b'a + Xa k can be reduced to a Hermitian form with any 
real A. Hence, in view of the assertion 1 we conclude that 



b' 2 t +b' 2 + (b' k + \) 2 >0, 



(29) 



S 



where A is an arbitrary real number. The above equality may be valid only when b' k = ^= G R. 
Choosing A = — b' k in (|2~9|) yields that b'f + 6' ■ > 0. Since 6'? + 6' 2 = (a x 6) 2 , hence we get 



3 — 1 J 

the desired inequality (ax 6) 2 > 0. The necessity is proved. 

In order to prove the sufficiency of the assertion 2, we consider transformation (123) with 



bidj — bjdi a k ab — b k a 2 



a 2 + a 2 v3? J of + a 2 



(30) 



This transformation leaves the matrix vo 2 ^ invariant, while the matrix b'a ( ffif ) transforms as 
follows 

r/,^ Vg x «6 
bo ->b a = = — o-j + -7=a k , (31) 

V a va 
whence it follows the sufficiency of the assertion 2. 

The proof of the assertion 3 of the lemma is similar to one of the assertion 2. The first three 



conditions are obtained with account of the assertion 2. A sequence of transformations (23) with 
a, b of the form (1271), ([SO]) transforms the matrix ca to become 



e ijk a(cxb) (axb)(axc) ac 
ca -»• c cr = ' = Oj + . Oj + ^=a fc . (32) 

\]{cxb) 2 \/(cxb) 2 Va~* Va 2 

Using the standard identities for the mixed vector products we establish that the coefficients by 
the matrices <jj, (jj, a k are real if and only if the relations 



jac, be, (a x b)c\ C R 



hold true. This completes the proof of Lemma 2. 

Lemma 2 plays the crucial role when reducing the potentials (PU[), to Hermitian forms. This 
is done as follows. Firstly, we reduce the QES operator to the Schrodinger form 

dy + f(y)aa + g(y)ba + h(y)ca + r(y). 

Here /, g, h, r are some linearly-independent scalar functions and a = (ai, a 2 , a 3 ), b = (6 1; b 2 , b 3 ), 
c= (ci, C2, C3) are complex constant vectors whose components depend on the parameters a, ft, 7. 
Next, using Lemma 2 we obtain the conditions for the parameters a, ft, 7 that provide a simulta- 
neous reducibility of the matrices aa, bo, ca to Hermitian forms. Then, making use of formulae 
(0), (0), (0) we find the form of the matrix A. Formulae Pfi|), ([31]) , (|52]) yield explicit forms 
of the transformed matrices aa, 6a, ca and, consequently, the Hermitian form of the matrix 
potential V(y). 



IV. QES matrix models 

Applying the algorithm mentioned at the end of the previous section we have obtained a complete 
description of QES matrix models ([17]) that can be reduced to Hermitian Schrodinger matrix 
operators. We list below the final result, namely, the restrictions on the choice of parameters 
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and the explicit forms of the QES Hermitian Schrodinger operators and then consider in some 
detail a derivation of the corresponding formulae for one of the six inequivalent cases. In the 
formulae below we denote the disjunction of two statements A and B as [A] \J[B]. 

Case 1. 71 = 72 = 73 = and 

[/3o, Pi, P2 e R] \J[P 2 = 0, pi = 2a 2 , p = a ± + i/i, /ieRj; 



H[y] = 



1 



{-Plx 4 - \2p x p 2 + Aa 2 p 2 (m - l)]x d 



+ 



4(a 2 x 2 + otix + ao) 
2a 2 p x - 2a x p 2 - p\ - 2p p 2 - A ai p 2 (m - 1) 



+ [Aa 2 p - 2p p 1 - Ama p 2 ] x + 2a!p - 2a p! - pi 
+Ap 2 (a 2 x 2 + aix + a )ai + (4/3 2 x + 25) (a 2 x 2 + a\x + a )a 3 } 
^a 2 3(2a 2 x + a±) 2 

2 16(a 2 x 2 + a\x + a ) 



x=f- 1 (y) 



A 



I. 



Case 2. p 2) 5 = and 

2a 2 p! - p\ G R, 2a 2 p - PoPi e R, 2a!/3 - 2/3i« - /?o - 7 2 e R, 

"(2a 2 - A) 2 7- > 0] V [2« 2 - A = 0] , [(c*i - A,) 2 7 2 > o] V - A) = 0] ; 



H[v\ 



P l {2a 2 -p l )x 2 + 2p {2a 2 -p 1 )x 



y ' \4(a 2 x 2 + a\x + ao) ^ 
+2 ai p - 2p iao - Pl - 7 2 + [2(2a 2 - p x )x + 2{a x - P )] ^fa 3 
a 2 3(2o; 2 x + aif 



+ 



2 16(a 2 x 2 + a±x + ccq) 



A 



Ai 2 (7i, 72)A 23 ( v / 7i +7l> 73), 7i + ll ^ 0. 



(If 7 2 + 75 = 0, then one can choose another matrix A (27) with 7 2 + 7 2 7^ 0). 
Case 3. a 2 ^ 0,p 2 ^ and 

{/3 2 ,7i} C Re, 73 = 0, 7 2 = y^ 2 - 2a 2 7i» «27i < 0, 

/3i = 2a 2 + /3 2 — , /3 = «i + /3 2 — 
« 2 a 2 



3(2a 2 a; + ai)" 



+ 



2_4 



2 16(a 2 :r 2 + ckix + «o) 4(a 2 :r 2 + ctix + a ) 
& 2 



2/3 2 2 — + 4a 2/ 3 2 m 

a 2 



X 



Q: 



(af + 2a «2) + 2ai/5 2 (l + 2m) 



x 
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2a 1 [3 2 (a 1 a 2 + a /3 2 ) 



a\ 



+ Aao(3 2 m 



2 r>2 a Cl An 

x + a, - p 2 ^ - 4/? 2 4:a a 2 - 2a 2 7i 

a 2 a 2 



sin y6(y)^-2a 2 r )i ) (J\ + cos ( 9(y)^-2a 2 ^i ) a 3 
(9{y)^2aVf 



+A(3 2 x(a 2 x 2 + d\x + ao) 

"sin 

+2(a 2 a^ + ot\x + ao) 



+ cos ( 6(y) x /-2a 2 j 1 ) ( J\ a2 =^~ a ^ + Sa 3 



V-2a 2 7i 



— ^y / -2a 2 7i°"i - 2/? 2 \/7i - 2a 2 7i°"3^ 



VV 7i V 7i 
Case 4. a 2 ^ 0, & = 0. 

Subcase 4.1. 5^0, 71,72 do not vanish simultaneously and 

7i - 72 < 0) 73 = ifJ>, {a 4 , S} C R, i(ai - f3 ) G R, /3i = 2a 2 ; 



j a 2 3(2a 2 ^ + «i) 2 



— { -fil + 2a 1/ 3 - 2a (3i - 7' 



2 16(a2# 2 + ®ix + a ) 4£ 

sin (e(y)^f^ 



+2(0:2^ + a±x + ao) 



, -^73 V72 - 7i*2 + ^(7i 2 - 7 2 2 )^3 / r— 2 
+ cos I 9{y)yJ--ft 
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+ 



25 'a-iD ' 3 



7 2 



2<fai7 3 (2e*i - 2 / 3 )7 J 2 + 2<5a 73 



-x + 



i\fl2 -ll°2 +730"3 



x=f~ 1 (y) 



11 li 
A = A 2 i (271,72). 
Subcase 4.2. 5^0, 71 = 72 = 0, 73 7^ and 
{5, P ± (2a 2 - /3i), #,(2a 2 - A), -/9 2 + 2c*iA> - 2a /3i, 7 3 (2a 2 - ft), 73(^1 - A))} C R; 



a 2 3(2o; 2 x + a\Y 



+ 



+ 

A = 1, 
Case 5. a 2 = 0, /3 2 7^ and 

,2 2 



2 16(o;2^ 2 + Cxi\x + ao) 4(a2a; 2 + a\x + ao) 
+2p (2a 2 - fc)x - Pi + 2«i/9o - 2(3^0 - 7 3 2 

25a 2 x 2 + 2rr((2a 2 - Pi)ls + Sa ± ) + 2(«i - /3 )73 + 25cd 03} 



{A(2a 2 -A)a; 2 



x=f~Hv) 



«i ^ 0, 71" - 72 < °> 7 2 < °> 73 



7 2 



2ai' 



{A), A, ft, 72, 5(7 2 2 - 7?) + 2^7173} C R, 

{i(2a /Vy3 - Pili + 2/3 2 ai7i + ^173), - A))7 2 + 2/3 2 a 7i + ^"073)} C R; 
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H[y] = 



9y + 



3aj 



+ 



(2 - 4m)a 1/ 3 2 - 0* - 2/5 /5 2 l a; 2 - [2/3 A + 4ma /3 2 ] a; 



+2ai / 9 - 2a A - /3 _ li + 4x(aia; + a ) 



sin (o(y)yf^H 



H rr; a 3 COS 

7i 

<K7 2 2 - 7i 2 ) + 2/3 27 i73 



•7 2 



+ 2(aix + a ) 



x 



+ 



+ 



-0"! 



2/3 2 727 2 



sin 



2&72 
72 - 7i 



=(7i + 



5(7 2 -7 2 2 )- 2/5 27 i73 



7 2 



-7/ 



(7i - iDlf 



ct 3 I cos ( 6{y)\J 7 2 



x- 



4a /j 2 7 3 - 2/317 2 + 4ai/3 2 7i + 2foi7 3 
if 



(2ai - 2/j )7 2 + 4a 0/ g 27l + 2fo 7 3 

~2 

A = A 2 i(i7i,7 2 )A 23 (^-^73^72 -7i,7 2 -72) 



x=f~ 1 (y) 



Case 6. a 2 = 0, /3 2 = 0. 

Subcase 6.1. 5^0, 71,72 do not vanish simultaneously and 

7 2 <0, {5 2 (7i- 7 2 2 ) <0, A), A}CR, 
{»(-A7. ? + <fai7s), »((«! - A))t? + ^073)} C R; 



A 



3« 2 



+ 



16(aia; + a ) 4(aix + ckq) 



- 2/5 0/ 9iX + 2<*iA> 



-2a /3i - /3q - 7 4 2 + 2(aix + a ) 

Kil - it) 



5^ 



ihi 



sin(%W-7 2 



^3 cos 9{y)J- ! y. 



\j {il - iDil 

(2ai - 2/3 )7 2 + 2^g 73 
if 



Z2 



+ 



X- 



\Rf 
2(3^1 + 2fo l73 

7 2 

J 

x=f~ 1 (y) 



= A 2 i(i7i, 72)A 23 ^73^72-71, 7i - t|) • 
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Subcase 6.2. 



7i = 72 = 0, 73 + 0, {Pi AA} C R, 

{-/3i73 + <fai, ("l - Po)l3 + 8a , -01 + 2a 1 p - 2a /?i} C R; 



H\ 



-\ 



3af 



1 16(«ix + «o 
1 



{-f3fx 2 - 2(3 Q (3 lX + 2ai/3 - 2a A - & - 7 3 



A 



4(«ix + «o) 

+2(«ix + a )[2x/3i(ai - 73) + 2(a x - /3 )j 3 + 2(3 1 a ](J?,} 
1. 



*=f~Hv) 



In the above formulae we denote the inverse of the function 



y = fix) 



^a 2 x 2 + a.\x + a ' 



(33) 



as / (y) and, what is more, the function 9 = 9{y) is defined as follows 



%) = - {/ 



dx 



a 2 ar + aiix + o>o 



(34) 



and 7 4 2 stands for 7? + 7| + 7|- 

The whole procedure of derivation of the above formulae is very cumbersome. That is why 
we restrict ourselves to indicating the principal steps of the derivation of the corresponding 
formulae for the case when a 2 7^ 0, /3 2 7^ omitting the secondary details. It is not difficult to 
prove that 7? 7^ 0. Indeed, suppose that the relation 7? = holds and consider the expression 
Q = U~ l {x)a^U{x) from fl2"U|). Making use of the Campbell-Hausdorff formula we get 

q2 

tt = 0-3 + 0(i7i<7 2 + 7 2 0i) - —l3li<?i-, 

where 9 is the function (0). Considering the coefficient at 9 2 , yields that 73 = (otherwise 
using Lemma 2 we get the inequality 7I7 2 7^ that contradicts to the assumption 7? = 0). Since 
the matrix coefficient at 9 has to be Hermitian, we get 7? = 7^ — 7! < 0. This contradiction 
proves that 7? 7^ 0. Taking into account the proved fact we represent the matrix potential ( |20| ) 
as follows 



v(v) 



a 2 

Y 



3(2ct 2 2 + oil) 2 
16(a 2 x 2 + ol\x + azo) 
1 



4(a; 2 ar + a\x + ao 



■A 



-1 



/3 2 V - [2/3^ 2 + 4a 2 /? 2 (m - l)]x d 



+ 



2a 2 pi - 2 ai p 2 - Pi - 2p p 2 - iatfoim - 1) 



x 
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+ [4a 2 /?o - 2/?o/?i - 4:ma (3 2 ] x + 2ai/3 - 2a /?i ~ Po ~ 7 



~,2 



-\-Ax(o>2X 2 + «iX + «o) 



folzill) 1 licr i + p 2 (l2(J-i + 271 02) {it) 1/2 sinh ( eJ^f 



+ [/52(-7i730"i - «72730"2 + (7i - tD^s)]^) 1 cosh (0\Jil ^ 



(35) 



(£72*1 + i(*7i - 2/3 2 7 3 )^2 - 2/3 27 2^)(7 2 )- 1/2 sinh 



+2(a;2X 2 + «ix + a ) 
+ [(2/3 2 (7 2 - 7 2 ) - £7173)0-1 - «(2/? 2 7i72 + 5 7273 )o- 2 + (5( 7 2 - 7 2 ) - 2/? 27l 7 3 )c73] 



(7 2 )" 1 cosh ( eJ^ 



+ \{-2(3 2 r t + 4a 2( 5 27 i + 25a 2 7 3 )x 2 + ((4a a - 2/3i)7 



^2 



+4ai/? 2 7i + 25ai 73 )x + 



(2ai - 2/3 ) 7 2 + 4a /3 27 i + 25a 7 3l/ ~ 



7 2 



](7 2 )- 1 W A 



where 6 1 = 0(y) is given by (|34|). 

Let us first suppose that 71,72 do not vanish simultaneously. We will prove that hence it 

follows that 7 2 G R. Consider the (non-zero) matrix coefficient at 4x£ cosh (^Q\J~^ ^ in the 
expression (|35| ) and suppose that = a + ib, with some non-zero real numbers a and 6. 
Now it is easy to prove that cosh [OyJ^ ^ = f(x) + ig(x), where f,g are linearly-independent 

real-valued functions. Considering the matrix coefficients of f(x), g(x) we see that in order to 
reduce the matrix ([35|) to a Hermitian form we should reduce to Hermitian forms the matrices 
A, iA which is impossible. This contradiction proves that 7? G R. 

sinh (0\J~^j; \ 



Consider next the non-zero matrix coefficients of Ax^ 
These coefficients can be represented in the form aa, ba, where 

a = #2(72, *7i>0)» b = /? 2 (-7i73, -^7273,7? - 7 2 ), 

and, what is more, 



4a£ cosh [eJjf in (M 



Applying Lemma 2 yields 



a x b = /3 2 (7! - 72 )(z 7l , -72,^73). 



& G R, 73 = i/i, /iGR, 7i - 7 2 < 0. 



Next we turn to the matrix coefficient of 2£ 



sinh (0\pjf 



which is of the form ca with 



c = {5j2-,i{$l\ — 2^73), — 2/3 2 7 2 ). Making use of the assertion 3 of Lemma 2 we obtain the 
conditions 

{7i, 72} CR, [ 7l = 0]\/[73 = 0]. 



14 



Considering in a similar way the matrix coefficient of 2£ cosh ( 9\/jf ) yields the following 
restrictions on the coefficients a, (3, 7: 



{/?2,7i} C R, 73 = 0,72 = V7i -2a 2 7i, "27i <0,Pi = 2a 2 + (3 2 — , 0o = "i + P2 — 

v a 2 a 2 

As a result we get the formulae of Case 2. 

One can prove in an analogous way that, provided 71 = 72 = 0,73 7^ 0, the matrix (]3 
cannot be reduced to a Hermitian form. 

V. Some examples. 

In this section we give several examples of Hermitian QES matrix Schrodinger operators that have 
a comparatively simple form and, furthermore, are in direct analogy to QES scalar Schrodinger 
operators. 

Example 1. Let us consider Case I of the previous section with ato = (3 2 — 1, the remaining 
coefficients being equal to zero. This choice of parameters yields the following Hermitian QES 
matrix Schrodinger operator: 

H[y] = dy - ^y 4 - my + a 3 y + o x . (36) 

The invariant space X of the above Schrodinger operator has the dimension 2m and is spanned 
by the vectors 

fj = exp (^j e x y\ g k = exp (^j (me 2 y k - ke y x k ~ l ), 

where j — 0, . . . , m — 2, k — 0, . . . , m, e\ = (1, 0) T , e 2 = (0, 1) T and m is an arbitrary natural 
number. 

Note that the basis vectors of the invariant space I are square-integrable on an interval 
(—00, B] with an arbitrary B < +00. It is also worth noting that there exists a QES scalar 
model of the same structure that has analogous properties ||. 

By construction, QES operator (|36| ) when restricted to the invariant space X becomes complex 
2m x 2m matrix M. However, the fact that operator (j36| ) is Hermitian does not guarantee that 
the matrix M will be Hermitian. It is straightforward to check that the necessary and sufficient 
conditions of hermiticity of the matrix M read as 

• basis vectors fj(y),gk(y) are square integrable on the interval [A, 5], 

• the condition 

(9 y rAy)yk(y) - ?Ay){dyn{y))\ B A = 0, (37) 

where r*j = /j, % = 0, . . . , m — 2 and r, = <fe_ m+ i, i = m — 1, . . . , 2m — 1, holds Vj, k = 
0, ...,2m- 1. 
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In the case considered relation ( p7|) does not hold and, consequently, the matrix M is not 
Hermitian. The next two examples are free of this drawback, since the basis vectors of their 
invariant spaces are square integrable on the interval (—00, +00). 

Example 2. Let us now consider Case I of the previous section with ai = l,/3 2 = — 1,/3 = 1/2, 
the remaining coefficients being equal to zero. This choice yields the following QES matrix 
Schrodinger operator 



H[y] 



y 256 



+ 



Am 



-y 



;y 2 <?3 



at. 



16 a 4' 

The invariant space X of this operator has the dimension 2m and is spanned by the vectors 



fi 



fjk 



exp 



exp 



-r 

64 



64 



2.7 



ei, 



m 'l) 



2k 



e 2 



2fc-2 



where j = 0, . . . , m — 2, k = 0, . . . , m. 

It is not difficult to verify that the basis vectors of the invariant space X are square integrable 
on the interval (—00, +00) and that the corresponding matrix M is Hermitian. One more 
remark is that there exists an analogous QES scalar Schrodinger operator whose invariant space 
has square integrable basis vectors (see, for more details 0, 

Example 3. Let us now consider Case III of the previous section with 0:2 = 1, /?2 = ~~ l>7i = — 1- 
This choice of parameters yields the following QES matrix Schrodinger operator: 



H[y] 



d y~\~ ^ ex P( -2 ^) + me w(-y) + ^exp(2y) 



+ 



+ 



rn ~~W~ sin(v / 2e s/ ) — cosf^v^e^) — exp(— y) sin(v / 2e y ) 
m — - — cos(v / 2e y ) H — — sin(-\/ / 2e ?/ ) — exp(— y) cos(v / 2e y ) 



0"! 



Furthermore, the invariant space X of this operator has the dimension 2m and is spanned by the 
vectors 



fi 
9k 



U 1 (y)exp(-jy)e u 

U' 1 (y) (mexp(-ky)e 2 - kexp(-(k - l)y)ei) 



where j = 0, 



m 



2, k = 0, . . . , m, m is an arbitrary natural number and 



U-\y) 



1 



: exp 



2^2 
x( v / 3 + V2-a 3 



exp 



cos(v / 2e i/ ) + 



iv3i 



cr 2 - 0i 



v/2 



sin(v / 2e y ) 



The basis vectors of the invariant space X are square integrable and the condition (|3 



holds. Indeed, the functions fj(y) and Qk{y) behave asymptotically as exp (- 



and 



16 



exp (— ^ 2fc ^ y ), correspondingly, with y — >■ +00. Furthermore, they behave as exp f- , 



(2j+l)y ^ 

x exp (— |e~ y )and exp ^_ ( 2fc + 1 )y ^ eX p ^_ |e~ y ), correspondingly, with y — > —00. This means 
that they vanish rapidly provided y — > ±00. 

Example 4- The last example to be presented here is the QES matrix model having a potential 
containing the Weierstrass function. To this end we consider the whole set of operators ([16]) and 
compose the Hamiltonian 

H[x] = D 2 + A 1 + 2B 2 . 
Reducing H[x] to the Schrodinger form yields the following QES matrix model: 



■run- u? l)w(y) 3 ^^ 2 - 1 ) 2 

16(w(y) 3 + w(y)) 



H[y] = d y + (to — to — l)w(y) 
2m 



w(yf 



l —{2cr 1 + {w{yf^?>w{y)) ( T Z 



Here to is an arbitrary natural number and w(y) is the Weierstrass function defined by the 
quadrature 

w(y) 



dx 



Vx 3 + 



X 



The invariant space X of the operator H[y] has the dimension 2m and is spanned by the 
vectors 

fj = (w(y) 3 + w(y)y^(l -ia 2 w(y))exp(-jy)e 1 , 
g k = (w(y) 3 + w{y)Y~i{l - ia 2 w(y)) (mexp(-ky)e 2 
-kexp(-(k - l)y)e x ) , 

where j = 0, . . . , m — 2, k = 0, . . . , to. 

Note that the first example of scalar QES model with elliptic potential has been constructed 
by Ushveridze (see || and the references therein). 



VI. Some conclusions 

A principal aim of the paper is to give a systematic algebraic treatment of Hermitian QES 
Hamiltonians within the framework of the approach to constructing QES matrix models sug- 
gested in our papers |L], [|. The whole procedure is based on a specific representation of the 
algebra o(2, 2) given by formulae @, (H), ©■ Making use of the fact that the algebra (0) has an 
infinite-dimensional invariant subspace (|5J) we have constructed in a systematic way six multi- 
parameter families of Hermitian QES Hamiltonians on line. Due to computational reasons we 
do not present here a systematic description of Hermitian QES Hamiltonians with potentials 
depending on elliptic functions (we give only an example of such Hamiltonian in Section V). 

The problem of constructing all Hermitian QES Hamiltonians of the form (|17]) having square- 
integrable eigenfunctions is also beyond the scope of the present paper. We restricted our analysis 
of this problem to giving two examples of such Hamiltonians postponing its further investigation 
for our future publications. 
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A very interesting problem is a comparison of the results of the present paper based on 
structure of representation space of the representation (f|), (|), (0) of the Lie algebra o(2, 2) to 



those of the paper [15], where some superalgebras of matrix-differential operators come into play. 



The link to the results of [JE5J is provided by the fact that the algebra o(2, 2) has a structure of 
a superalgebra. This is a consequence of the fact that operators (0) fulfill identities (Q). 

One more challenging problem is a utilization of the obtained results for integrating multi- 
dimensional Pauli equation with the help of the method of separation of variables. As an 
intermediate problem to be solved within the framework of the method in question is a reduction 
of the Pauli equation to four second-order systems of ordinary differential equations with the help 
of a separation Ansatz. The next step is studying whether the corresponding matrix-differential 
operators belong to one of the six classes of QES Hamiltonians constructed in Section IV. 

Investigation of the above enumerated problems is in progress now and we hope to report 
the results obtained in one of our future publications. 
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